Mechanisms of nonlinear spin-wave emission from a microwave driven nanocontact 
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We present a micromagnetic study of linear and nonlinear spin-wave modes excited in an extended 
permalloy thin film by a microwave driven nanocontact. We show that the linear mode having 
the frequency equal to the excitation frequency (/) is driven by the ac Oersted field component 
perpendicular to the static external field (applied in-plane of the sample) . The nonlinear mode with 
the frequency //2 is excited as an independent eigenmode within a parametric longitudinal pumping 
process (due ac Oersted field component parallel to the bias field). Spectral positions of those modes 
are determined both in the space and phase domain. The results are important for the transfer of 
information coded into spin-waves between nanocontacts, and for synchronization of spin transfer 
torque nano-oscillators. 

PACS numbers: 75.30.Ds, 75.78.Cd, 75.75.+a 



Spin waves in meso- and micro-sized magnetic struc- 
tures have attracted growing attention due to the 
prospective applications in magnon spintronic devices. 
For example, the use of phase and amplitude of spin 
waves as an information carrier in magnetic logic cir- 
cuits is extensively studied.— ~— The problem of spin- wave 
excitation is of crucial importance for the miniaturiza- 
tion of such devices to the nano-scale. The advanced 
research on device miniaturization makes it possible to 
study the spin waves exited by microwave driven nano- 
sized antennas (nanocontacts) £~— This task is not trivial: 
if the density of the microwave current is high enough, the 
spin- wave system is pushed far beyond the linear regime, 
and nonlinear processes appear. Moreover, in nano-sized 
structures, the properties like the threshold power and 
frequency limit of the nonlinear excited spin-wave modes 
can be tuned by a superimposed direct current flowing 
through the contact^ Vice versa, the influence of the mi- 
crowave currents on the magnetic oscillations caused by 
the spin transfer torque effect (STT) in direct current 
driven nanocontacts was demonstrated^*^ It was shown 
that the propagation length of the STT emitted spin 
waves can be enhanced by using an additional microwave 
current injected through the contact.— The microwave 
current plays also an important role in the synchroniza- 
tion of STT based nano-oscillators (STNO)^iS^ All 
these tasks need a deeper understanding of the char- 
acteristics of spin waves emitted by microwave driven 
nanocontacts. 

In this work we performed systematic micromagnetic 
simulations^ on magnetization dynamics in a nanocon- 
tact structure to reveal the underlaying mechanism of 
magnons excitation and to characterize the linear and 
nonlinear spin-wave modes detected experimentally in 
such devices 4 

In the simulations we considered a thin permalloy layer 
of 7 nm thickness with lateral dimensions of 1.5 x 1.5 /^m 2 



(see Fig.QJa)). As in the experiment^, an in-plane mag- 
netic field of 19.89 kA/m (250 Oe) was applied along 
the x-axis. A microwave current is injected through 
the nanocontact placed in the center of the structure 
as defined in Fig. [Ha). The resulting alternating Oer- 
sted field excites magnetization precession in the sys- 
tem. The circular shape contact has a diameter of 
100 nm. The standard material parameters of permal- 
loy used to simulate the magnetization dynamics are: 
saturation magnetization HqMs = 1 T, exchange stiff- 
ness constant A = 1.3 x 10~ n J/m and zero magneto- 
cystalline anisotropy was assumed. The simulated area 
was discretized into N x X N y X N z = 300 x 300 x 2 
cells, each cell having a size of 5 x 5 x 3.5 nm 3 . Ini- 
tially, the magnetization is aligned in-plane along the 
x-axis. To model the parametric processes in our sys- 
tem, we considered a random fluctuation of the magne- 
tization within 0.5° with respect to the x-axis, which is 
caused by thermal noise in real systems. The damping 
boundary conditions were used to avoid spin-wave re- 
flection at the edges of the simulation area. A spatial 
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FIG. 1. (a) Scheme of the sample with the contact in the 
center and a nonmagnetic metallic asymmetric lead situated 
at 60 nm above the magnetic layer, (b) Spatial distribution 
of the Oersted field components parallel and perpendicular 
(inset) to the external applied field. The inset figure has a 
lateral size of 300 nm x 300 nm. 
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distribution of the damping parameter (a(x, y)) in the 
Landau-Lifshitz-Gilbert equation (with the STT term in- 
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eluded) given by: a = ao + A 

with A = 1, R = 800 nm, and a = 100 nm ensures a 
constant value ao = 0.01 in the center of the sample for 
a radius of ~ 500 nm. The variables x and y denote 
the spatial coordinates. With the above parameters, the 
damping is increased more than fifty times at the edges 
of the simulation area, thus the spin-wave intensity de- 
cays strongly towards the boundaries. We calculate^ 
the current distribution and the corresponding Oersted 
field in the Permalloy layer with and without the contri- 
bution of an asymmetric top electrode (according to the 
real structure in the experiment^) . It was found that the 
lead induces an asymmetry for the Oersted field compo- 
nent parallel to the external field while the perpendicular 
component remains virtually unchanged (Fig. [lib)). I R 
the presented simulations we used the distribution of the 
Oersted field with the contribution of the lead included. 

From the experiment 8 it was found that two conditions 
must be fulfilled in order to generate a spin-wave mode 
with half of the excitation frequency. Firstly, the excita- 
tion frequency must be at least twice the frequency at the 
bottom of the spin- wave dispersion band (noted / ), and 
secondly, the applied microwave power is above a certain 
threshold value. In our configuration f Q ~ 4.55 GHz, 
we set the excitation frequency to 11 GHz to satisfy the 
first condition. The exact value for the threshold power 
could not be determined by simulations due to extremely 
long computation times. We used an applied microwave 
power sufficiently large in order to observe the nonlinear 
excited spin waves with good contrast. 

Analyzing the magnetization fluctuation in time for 
each computational cell and performing the Fast Fourier 
Transform (FFT), we obtained the frequency spectra for 
each component of the magnetization. The behavior of 
the magnetization in one computational cell (top edge of 
the contact) is presented in Fig. [5] One sees a response 
with the same frequency as the excitation one (/), but 
also the generation of nonlinear modes: 0.5/, 1.5/ and 2/ 
(Fig. H^a)). Furthermore, the precession of the magne- 
tization vector tip has a steady complicated trajectory 
when the nonlinear modes are generated (see Fig. [2Kb)), 
different from typical clamshell orbits found in thin films 
for linear excitations. An additional small loop was ob- 
served on the main gyration motion. This peculiarity 
originates from the combination of partial magnetization 
motions at different frequencies. 

Our main task is to identify the excitation mechanism, 
spatial localization and the wavevector directions as well 
as wave numbers for both linear and nonlinear excited 
spin- waves modes. 

First, we address the generation of the linear mode for 
which the frequency matches the microwave one (/). The 
spatial distribution is obtained by taking the amplitude 
value at / = 11 GHz from the frequency spectra of each 
computational cell and displaying them as a function of 



the position. The obtained map is shown in Fig. . 
The black circle in the center represents the position of 
the nanocontact. One can observe two maxima that are 
located on the left and right side edges of the contact. 
This can be understood by taking into account that the 
main driving field at the excitation frequency is the Oer- 
sted field component perpendicular to the bias field, with 
a very similar amplitude distribution in space (see the 
inset from Fig. [2(a)). From Fig. [3ja), we identified two 
types of waves excited by the nanocontact: travelling 
and non-propagating waves. The travelling spin waves 
exhibit in four beams starting from those two maxima 
regions, and are propagating along directions nearly per- 
pendicular to the bias field (applied on the x-axis), while 
the non-propagating ones can be observed in the rest of 
space. Moreover, the two maxima regions are acting as 
two independent sources which emit spin waves out of 
phase (180°) in respect to each other. This is illustrated 
on the snapshot image of the magnetization oscillation 
pattern in the entire permalloy layer (Fig. 03h)) with 
the M z component color coded. Since the two spin-wave 
beams that are travelling on the same side of the contact 
(up or down) are out of phase, a destructive interference 
takes place between them. Consequently, the spin-wave 
intensity is strongly reduced on the direction perpendicu- 
lar to the bias field in the center of the contact (Figs. Eta) 
and Mb))- 

The wavevectors of the emitted spin waves are deter- 
mined by applying the two-dimensional Fourier trans- 
form in space to the already calculated components (real 
and imaginary parts) of the temporal FFT for each com- 
putational cell. Figure M, c ) shows the distribution of 
the / mode in the wavevector space. k x and k y axis 
are the wave numbers that indicate the in-plane direc- 
tions of propagation parallel and perpendicular to the 
external applied magnetic field, respectively. The black 
and yellow contours in the figure represent the analytical 
constant-frequency curve calculated from the dispersion 
relation^ for a frequency of 11 GHz and 5.5 GHz, respec- 
tively. The wavevector data obtained from the simulation 
is in very good agreement with the analytical constant- 
frequency curve calculated at 11 GHz (the black line in 
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FIG. 2. (a) Magnetization oscillation spectrum for a cell 
at the top edge of the contact at the excitation frequency 
/ = 11 GHz. (b) Trajectory of the local magnetization vector 
tip for the same cell at the edge of the contact. 
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the figure). 

The analysis in the wavevector space (Fig. |3{c)) shows 
that indeed the most efficiently excited spin waves are the 
ones that are propagating with wavevectors nearly per- 
pendicular with respect to the bias field. The excitation 
efficiency is decreasing fast with the spin- wave propaga- 
tion angle, becoming close to zero for spin waves with 
wavevectors parallel to the bias field. This is associated 
with the known effect of decreasing excitation efficiency 
for spin- waves with short wavelengths. Nevertheless, the 
extracted data from the simulation shows another two 
strong maxima located at small wave numbers. Those 
maxima are related to the forced non-resonant driving 
of magnetization by the excitation field. Their location 
is traced back to the spatial configuration of the driv- 
ing ac field. Figure |3][d) presents the distribution of the 
perpendicular component of the Oersted field computed 
in the wavevector space. By comparing these distribu- 
tions (magnetization and excitation field) one observes 
that the maxima are located at the same wave numbers. 
Thus, one can see that the perpendicular component of 
the driving Oersted field is responsible for excitation of 
both propagating and localized spin-wave modes at the 
signal frequency /. The excitation efficiency for travelling 
waves reaches the maximum at directions perpendicular 
to the bias field. 

Next we discuss the process of excitation and the char- 
acterization of the non-linear spin- wave modes, especially 



E 

■E 



(a) >H 

m 



-400 -200 200 

x [nm] 



high 





(d) 



I high 



-7.0x10 7 -3.5x10' 0.0 3.5x10' 7.0x10 7 -7.0x10 7 -3.5x10' 0.0 3.5x10' 7.0x10' 

K [m ] k„ [m 1 ] 

FIG. 3. (a) Spatial distribution of the / mode and the position 
of the contact in the center (black circle), (b) Snapshot image 
of the magnetization oscillation pattern with the component 
M z color coded (red - positive values, blue - negative values). 

(c) Distribution of the / mode in the wavevector space. The 
black and yellow lines represent the constant-frequency curves 
calculated analytical for 11 GHz and 5.5 GHz, respectively. 

(d) Distribution of the perpendicular component of the Oer- 
sted field computed in the wavevector space. 



the f/2 mode, is described in the second part of the ar- 
ticle. 

The appearance of the f/2 mode (= 5.5 GHz), can be 
potentially attributed to three different excitation mech- 
anisms. Mechanism 1 is related to the three-magnon 
scattering process (confluence and splitting of magnons - 
quanta of the spin waves which are eigenexcitation of the 
magnetic system). Within mechanism 2, the H y com- 
ponent of the driving Oersted field excites an ac mag- 
netization perpendicular to the bias field. This ac mag- 
netization serves as a pump for the //2 mode, via the 
so-called process of parametric excitation of spin waves 
under perpendicular pumping. Within this process, the 
f/2 mode can be generated directly from the forced ex- 
citation of magnetization at the frequency /. Mechanism 
3 is a direct excitation of the f/2 mode by the Oersted 
field component parallel to the bias field, a process called 
longitudinal (or parallel) parametric pumping^ 

The radiation pattern extracted from the simulation 
for the f/2 mode is presented in Fig. 2{a). The black 
circle indicates the position of the nanocontact. One can 
observe that the f/2 mode is localized outside the con- 
tact area in the corresponding regions of travelling spin 
waves for the directly excited / mode. This observa- 
tion may suggest that we are dealing with a splitting of 
magnons of the / mode. The law of energy conservation 
for such processes is fulfilled since any mode could be ob- 
tained as a combination (confluence or splitting) of the 
other modes (see Fig. HJa)). However, it has to be noted 
that the wave numbers extracted from the magnetization 
intensity distributions in the phase space for both / and 
f/2 (Fig.G^c) and Fig.H») modes do not fulfill the mo- 
mentum conservation law. Moreover, from the constant- 
frequency curves calculated analytically for 11 GHz and 
5.5 GHz (black and yellow curves from Fig.^c)) one sees 
that the conservation of momentum during the splitting 
process is not satisfied for any combination of wavevec- 
tors. In other words, there is no pair of wavevectors at 
the frequency of f/2 — 5.5 GHz that can combine into 
one corresponding to 11 GHz, or vice versa. This strongly 
contradicts the assumption that three-magnon scattering 
processes take place in the given configuration, ruling out 
mechanism 1. The analysis in the phase space of the f/2 
mode (displayed in Fig.|^b)) shows a propagation of spin 
waves with wavevectors under a maximum angle of ^15° 
with respect to the bias field (x-direction) . The asym- 
metry is related to additional contributions of the lead 
to the internal field (similar for the asymmetry in the 
radiation pattern from Fig. [4ja)). 

In order to discriminate between mechanism 2 and 3 
for the f/2 mode generation (perpendicular and parallel 
parametric pumping), we performed two additional sim- 
ulations, in each one taking into account only one compo- 
nent of the Oersted field, H y and H x , respectively. The 
applied power was kept at the same value as previously. 

By applying only the perpendicular component H y , 
the qualitative behavior for the / mode remains virtually 
unchanged in both space and phase domain. Practically 
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FIG. 4. Spatial distribution of the f/2 mode obtained by tak- 
ing into account in the simulations: (a) both components (H x , 
H y ) of the Oersted field; (d) only the parallel component H x . 
The black circle indicates the position of the nanocontact. (b) 
The extracted distribution of the f/2 mode in the wavevector 
space considering components (H x , H y ) of the Oersted field; 
(c) Magnetization oscillation spectrum for the cell at the top 
edge of the contact considering i) only H x component (black 
line) and ii) only H y component (red line) distributions of the 
ac Oersted field. 



no difference was observed in comparison with the two 
distributions presented in Fig. [SJa) and Fig. (3Jc). This 
confirms that the / mode derives from the perpendicular 
component of the ac Oersted field. However, from the fre- 
quency spectra (red curve in Fig. Etc)) it was observed 
that the f/2 mode had vanished. Even if the pump- 
ing power was doubled, no peak was seen at the corre- 
sponding frequency for the f/2 mode. Therefore, one can 
conclude that perpendicular parametric pumping (mech- 
anism 2) is not responsible for the generation of the //2 
mode. Also, it evidences that the splitting of magnons 
does not occur even if there is a strong excitation of the 
/ mode. 

In the longitudinal parametric pumping process (mech- 
anism 3), the magnetization is excited by an ac field par- 



allel to the steady one. Within this process, spin waves 
having half of the excitation frequency are amplified di- 
rectly from the thermal bath. The frequency spectra ob- 
tained from the simulation (black line in Fig. Etc)) show 
two peaks corresponding to the / /2 and 3/ /2 modes and 
no intensity for the / mode. The spatial distribution of 
the //2 mode is presented in Fig.|4jd). Two maxima are 
observed and are consistent with the spatial distribution 
of the excitation field (only H x component). The strong 
asymmetry is derived from the additional field created 
by the lead on top of the sample. One can notice that 
the intensity distribution of the f/2 mode obtained in 
this case is very similar to the one resulting when both 
components of the Oersted field were taken into account 
(Fig. Eta)). Moreover, the analysis of the f/2 mode in 
the phase space shows exactly the same distribution of 
the wave numbers as was found when both H x and H y 
components were considered (Fig.Etb)). The slight shift 
observed by comparing the two magnetization distribu- 
tions from Fig. Eta) and Fig. Etc) can be interpreted as a 
distortion caused by the travelling waves of the / mode 
on the f/2 mode. Altogether, this shows clearly that 
the f/2 mode is excited as an independent eigenmode by 
parametric longitudinal pumping. 

In conclusion, we performed a micromagnetic study of 
the linear and nonlinear spin- wave modes excited by mi- 
crowave currents in a nanocontact structure. We have 
shown that the dominant spin-wave mode is driven by 
the dynamical Oersted field component perpendicular to 
the static external field. The nonlinear f/2 mode derives 
from the parametric longitudinal pumping effect (due al- 
ternating Oersted field component parallel to the bias 
field) . Spectral positions of those modes are determined 
both in the space and frequency domains. The / mode 
emitted from the points localized at edge of the con- 
tact is propagating nearly perpendicular to the bias field. 
The f/2 mode is localized outside the contact area and 
is propagating under a small angle with respect to the 
steady field. These results are important for the realiza- 
tion of nano-sized spin-wave antennas in magnon based 
spintronic devices and can be used for the design of spin- 
wave coupled arrays of nano-oscillators. 
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